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ABSTRACT 

Context. High-mass X-ray binaries (HMXBs) might have contributed a non-negligible fraction of the energy feedback to the interstel¬ 
lar and intergalactic media at high redshift, becoming important sources for the heating and ionization history of the Universe. How¬ 
ever, the importance of this contribution depends on the hypothesized increase in the number of HMXBs formed in low-metallicity 
galaxies and in their luminosities. 

Aims. In this work we test the aforementioned hypothesis, and quantify the metallicity dependence of HMXB population properties. 
Methods. We compile from the literature a large set of data on the sizes and X-ray luminosities of HMXB populations in nearby 
galaxies with known metallicities and star formation rates. We use Bayesian inference to fit simple Monte Carlo models that describe 
the metallicity dependence of the size and luminosity of the HMXB populations. 

Results. We find that HMXBs are typically ten times more numerous per unit star formation rate in low-metallicity galaxies (12 + 
log(0/H) < 8, namely < 20% solar) than in solar-metallicity galaxies. The metallicity dependence of the luminosity of HMXBs is 
small compared to that of the population size. 

Conclusions. Our results support the hypothesis that HMXBs are more numerous in low-metallicity galaxies, implying the need to 
investigate the feedback in the form of X-rays and energetic mass outflows of these high-energy sources during cosmic dawn. 
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1. Introduction 

It has been proposed that the energy feedback of high-mass X- 
ray binaries (HMXBs) plays a major role in the thermal and ion¬ 
ization histories of the Universe (Power et al.|2009, Mirabel et| 


al.|201 l[[Jeon et al.|2013[|Power et a 1, 1 201 3[ Knevitt et al. 2014| >, 

and that it significantly affects the cosmic star formation rate 


and the evolution of low-mass galaxies (Justham & Schawin- 
ski|2012||Artale et al.|2015| ). The importance of these sources is 
based on the hypothesis that they are more numerous, and possi¬ 
bly more luminous, in low-metallicity stellar populations. This, 
coupled with the chemical evolution of the Universe, suggests 
that in the early Universe the contribution of HMXBs to the en¬ 
ergy feedback into the interstellar and intergalactic media (ISM 
and IGM) would have been at least as large as that of other high- 
energy feedback sources, such as supernovae. 

However, direct evidence for the increase in the number and 
luminosity of HMXBs in low-metallicity stellar populations is 
still scarce. [Grimm et aT7| ( |2003| > and |Mineo et al. | ( |20 1 2[ ) have car¬ 
ried out exhaustive observational studies, compiling statistically 
significant samples of HMXBs in local galaxies. They found that 
the total X-ray luminosity of a HMXB population (Lx) scales 
with the star formation rate (SFR) of its host galaxy, which at 
first sight seems to imply that there are no metallicity effects. 
However, their sample includes mainly high-SFR galaxies that, 
owing to the correlation between SFR and metallicity, would be 
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high-metallicity systems. Hence, no metallicity effects are ex¬ 
pected. Trying to overcome this limitation, Kaaret et al.| (|2011 1 
and|Brorby et al~| ( |20 1 4~| studied a sample of local low-SFR, low- 
metallicity blue compact dwarfs (BCDs). They found that these 
galaxies host HMXB populations ten times more numerous per 
unit SFR, than galaxies in the sample of |Mineo et al.| ( |20T2| . 

On the other hand, Basu-Zych et al. ( 2013a] ) have shown that 
the Lx-SFR relation evolves with redshift, which could be due 
to metallicity effects. Fragos et al. ( |2013[ ) used population syn¬ 
thesis models coupled to large scale cosmological simulations 
to investigate the emission of HMXB populations at different 
redshifts. They found that the integrated X-ray luminosity per 
unit SFR increases with redshift, exceeding that of other X-ray 
sources (low-mass X-ray binaries, active galactic nuclei) at high 
redshifts. In addition, Basu-Zych et al. ( 201 3b[ > showed that the 
relation between the ratio Lx /SFR and the gas-phase metallicity 
of a sample of z < 0.1 Lyman-break analogs with slightly sub¬ 
solar metallicities (12 + log(0/H) > 8.1) is consistent with the 
theoretical expectations of Fragos et al. ( 2013| >. 

Although previous works showed the existence of a metal¬ 
licity dependence of HMXB populations, a discrimination of its 
effects on the size of these populations from those on the intrinsic 
luminosities is still lacking. To explore this issue, in this work we 
compile from the literature a large sample of local galaxies for 
which individual HMXBs can be resolved. Our sample extends 
those of previous works (Mi neo et al.|2012] Brorby et al. [ 2014) 
by including homogeneous estimations of the metallicities of the 
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galaxies, hence allowing the study of a larger metallicity range. 
Moreover, our analysis extends the work of previous authors by 
requiring our models to match both the data of the HMXB pop¬ 
ulation sizes and luminosities. In this paper we present the con¬ 
struction of the observational sample (Sect. |5J, and investigate 
the dependence on metallicity of the properties of HMXB popu¬ 
lations (Sect, 
sions in Sect. 


3l». We discuss the results and present our conclu- 

m 


2. The sample 

In order to investigate the dependence of the properties of 
HMXB populations on the metallicity of their progenitors, for 
49 galaxies in the Local Universe we compiled from the liter¬ 
ature the number of binaries, the integrated X-ray luminosities, 


star formation rates (SFR) and metallicities. Following Mineo 
et al. \2b\2), we selected galaxies that satisfy two criteria: a) 


high specific star formation rate (sSFR = SFR/M*), and b) small 
distance. As the production of HMXBs depends on the SFR of 


the galaxy (Grimm et al. 2003!, and the number of low-mass 
X-ray binaries is proportional to its stellar mass M„ ( jGilfanovj 
|2004[ ), the first criterion is useful to select galaxies in which 
HMXBs dominate the population of X-ray binaries. In partic¬ 
ular, we chose sSFR > 10 10 yr '. The distance criterion was 
imposed because we are interested in galaxies where individ¬ 
ual HMXBs can be resolved. The galaxies in our sample were 
chosen by requiring that their distance is D < 65 Mpc, or their 
redshift is z < 0.015. 

The main source of objects for our sample is the work of 
|Mineo et al.| ( |2012) >. This work compiles the number of sources, 
luminosities (both of individual XRBs and the total integrated 
luminosity of the galaxy), and star formation rates for 29 nearby 
(D < 40Mpc), star-forming galaxies with SFR ~ 0.1 - 20 M e . 
The sample of these authors is based on a homogeneous set of 
observations in X-rays, infrared and UV from Chandra , Spitzer, 
GALEX and 2MASS, with very specific selection criteria. We 
extended their sample in two ways. First, we compiled from the 
literature the metallicities of the galaxies when available. This 
left us with only 19 galaxies, most of them being high-metallicity 
systems. As a second step we added to this sample a set of low- 
metallicity galaxies because the investigation of a large abun¬ 
dance range is needed to assess clearly the effects of metallic¬ 
ity on HMXB populations. The low-metallicity galaxies come 
mainly from the work of |Brorby et al. | (20141, who studied a 
set of 25 nearby (D < 30 Mpc) low-metallicity BCDs. Another 
5 galaxies from differen t works (|Ghigo et al.|1983| |Thuan et al.| 
2004; |Winter et al.|2006||Kaaret et al.|201 1[ > were also added. We 

note that 20 of the 49 galaxies in our sample have no HMXBs 
detected; however these data still provide useful constraints on 
the metallicity dependence, as we will show in Appendix[A] The 
compilation of abundances, together with the extension of the 
metallicity range and the use of both the observed population 
sizes and integrated luminosities, would allow us to go beyond 
the results of previous works. As the low- and high-metallicity 
samples were taken from different sources, with different un¬ 
certainties and selection effects, a standardization of the data 
is needed to make statistical comparisons and draw meaningful 
conclusions. In the following sections we describe the steps fol¬ 
lowed to make the data homogeneous, and the caveats involved 
in the process. 
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2.1. Number of binaries and their luminosity 

As all the works from which we construct our sample give the 
data of individual sources in the galaxy, we took directly the 
observed number of sources N and their individual X-ray lu¬ 
minosities from the aforementioned works (the complete data 
of our sample is given in Table [TJ. The luminosity threshold 
L* above which sources are detected was also recorded. Al¬ 
though we cannot ensure that all the sources listed are indeed 
HMXBs, the selection criteria imposed minimize the contami¬ 
nation by other types of sources. In the case of the data taken 
from Mineo et al. ( |2012j >, we only used objects with location 
flag equal to 1 in their Table Al, which are those sources ly¬ 
ing in the HMXB-dominated regions of the galaxies. For the 
data of low-metallicity galaxies, the association of the few X-ray 
sources found with star-forming regions makes highly improba¬ 
ble of them being other kind of sources. 

X-ray luminosities of individual sources are given in differ¬ 
ent energy ranges in the literature, according to the instrument 
used to observe them. To make the data homogeneous, we cor¬ 
rected all the observations to the 0.5 - 8 keV band, which is that 
of most observations of our sample. To perform the correction, 
we assumed that the X-ray spectrum of individual sources is a 
power-law, and used for each source the spectral index reported 
in the literature. In the case it is not available, we took the stan¬ 
dard value —1.7 ( Fabbiano{2006| l. The luminosity thresholds of 
the observations were corrected accordingly. Throughout this pa¬ 
per, we will express all the luminosities in the 0.5 - 8 keV band. 

The total X-ray luminosity L ( x ' of each galaxy in our sample 
is taken as the sum of the luminosities L, of individual sources 
above L t j,. However, the latter differs for each galaxy. Aiming at 
obtaining a homogeneous luminosity sample and following [Mr] 
neo et al](j2012|, we also computed the total luminosities of the 
galaxies corrected to a common threshold L t h,i = 10 36 erg s '.In 
the case of NGC 1569, as L t h < L t h,i- this correction was made 
by subtracting the luminosity of the sources below L t h,i- For the 
remaining galaxies, we used a power-law X-ray luminosity func¬ 
tion (XLF) with the index (a = -1.58) and normalization factor 


(q = 1.49 (10 38 erg s I ) Op58 M 0 I yr) given by Mineo et al. (2012 1 , 
hence obtaining the corrected luminosity 


r G,36 _ t G 


X 


= L% + q SFR 


f*. 

‘-'^th,! 


r (T+1 


dL. 


(1) 


In most cases, the integral in the previous equation amounts to 
less than 11% of the final value, except for UGC 5720 (19%), 
NGC 4194 (28%), and NGC 3079 (29%). We conservatively as¬ 
sumed here that the XLF does not depend on metallicity. This 
choice is not crucial because the corrections are small, owing to 
the weak dependence of the integrated luminosity on L th . How¬ 


ever, if the higher normalization factor given by Brorby et al. 

' ( |2014[ > were used for low-metallicity galaxies, a slightly larger 
total luminosity would have been found, making our conclusions 
stronger. 

It should be pointed out that the aforementioned correction 
on the luminosities was not applied to the galaxies without de¬ 
tected sources. In this case, any trend found in the corrected val¬ 
ues would reflect the corresponding trends of the model used in 
Eqn.[T|for the XLF, rather than those of the original data, hence 
biasing the results. 

Following the above reasoning, we also corrected the num¬ 
ber of sources detected in each galaxy to a common threshold 
Eth ,2 = 10 38 ergs _1 . The choice of this threshold value was dic¬ 
tated by the fact that the observed number of sources does de- 
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Galaxy 

D 

[Mpc] 

12 + log(0/H) 

SFR 

[Moyr 1 ] 

log 

[erg s~ ] 

log L * 
[ergs^ 1 ] 

N 

1 T G,36 

log 

[erg s- 1 ] 

1 rG,38 

log L x 
[ergs- 1 ] 
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2 
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- 

0 
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1 
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1 
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39.55 

2 
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1 
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1 
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38.45 

1 

39.00 
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1 
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20.8 

7.51 

0.040 
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38.52 

1 
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1 
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8.67 

7.64 

0.015 

38.34 

38.11 

1 

38.35 

38.34 

1 

RC2 A 1228+12'' 

21.2 

7.64 

0.031 

38.53 

38.53 

1 

38.55 

38.54 

1 
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0.014 

38.62 

37.14 

2 
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1 
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1 
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23 
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6 
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35.71 

10 
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- 

0 
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0.44 
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8 
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4 
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17 
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1 
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4 
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- 

0 
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36.16 

42 
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2 
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38.41 

36.19 

14 

38.42 

38.17 

1 

NGC 4490" 

7.8 

8.36' 

1.8 

40.06 

37.1 

32 

40.07 

40.03 

11 

NGC 3034" 

3.9 

8.36“ 

10.5 

39.86 

36.86 

54 

39.90 

39.75 

12 

NGC 4038-39" 

13.8 

8.4' 

5.4 

40.25 

36.92 

83 

40.26 

40.20 

20 

NGC 7793" 

4 

00 

0.29 

38.36 

36.55 

9 

38.38 

38.20 

1 

NGC 3310" 

19.8 

8.44' 

7.1 
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37.84 

23 

40.79 

40.77 

22 
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24.9 
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39.57 

38.36 

4 

39.66 

39.60 

6 

NGC 5457" 

6.7 

8.55" 

1.5 

39.60 

36.36 

96 

39.61 

39.49 

11 

NGC 3079" 

18.2 

8.57 ; ' 

6 

39.65 

37.98 

14 

39.80 

39.65 

14 

NGC 5194" 

7.6 
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40.01 
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69 

40.02 

39.91 

16 
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25 
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40.52 

25 
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39.1 
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38.64 

4 
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40.49 

29 
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12 
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37.74 

25 

40.49 

40.45 

22 
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7.16 
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38.21 

0 

- 
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0 
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7.24 

0.012 

- 

38.38 

0 

- 

- 

0 
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7.26 
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- 
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0 
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24.6 
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38.53 
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- 

0 

6 dFJ 0405204-364859'' 

11 

7.34 

0.013 

- 

38.33 

0 

- 

- 

0 

HS 0822+3542' 

12.7 

7.35 

0.004 

- 

38.48 

0 

- 

- 

0 

J 085946.92+392305.6' 

10.9 

7.45 

0.002 

- 

38.36 

0 

- 

- 

0 

KUG 0937+298'' 

11.2 

7.45 

0.015 

- 

38.36 

0 

- 

- 

0 

UGC4483' 

3.44 

7.54 

0.004 

- 

37.55 

0 

- 

- 

0 

J 120122.3+021108.5' 

18.4 

7.55 

0.008 

- 

38.58 

0 

- 

- 

0 

KUG 0201-103'' 

22.7 

7.56 

0.015 

- 

38.54 

0 

- 

- 

0 

KUG 1013+381'' 

19.6 

7.58 

0.073 

- 

38.56 

0 

- 

- 

0 

SBS 1415+437'' 

13.7 

7.59 

0.040 

- 

38.52 

0 

- 

- 

0 

SBS 1102+606' 

19.9 

7.64 

0.035 

- 

38.54 

0 

- 

- 

0 

KUG 0942+551'' 

24.4 

7.66 
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- 

38.52 

0 

- 

- 

0 

KUG 0743+513'' 

8.6 

7.68 
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- 

38.14 

0 

- 

- 

0 
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- 

37.96 

0 

- 
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0 
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- 

37.66 

0 

- 

- 

0 
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8 * 

1.32”* 

- 
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0 

- 
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Notes. References: (i) |Brorby et al. (2014), (,/) Mineo et al. (|2012), ( '" ) Kaaret et al.|(|2011), (a) Izotov & Thuan (1999 

, (w |lzotov& Thuan (2011), 

Calzetti et al. (2007], ^Kobulnicky & Ski lman (1997), 1/1 Lopez-Sanchez & Esteban (2010), Esteban et al. 

2009 ), ( Engelbracht et al. 

(2008), ^Pilyugin & Thuan (2007 iT” Mirabel et al. (1992), ^ Bibby & Crowther (2010)7” Denicolo et al. (2002 

WT 

Hirashita et al. ( 2002 ), 

w Cedres et al. (2004),"^’ Bresolin et al. (2004),"^ Werk et al. (2011), w Storchi-Bergman et al. (1996 ), {r) Kniazev et al. 

(2004 ), 171 Thuan et al. 

('2004)7* Winter et al. (2006iT nrT lGhigo et al. ( 1983], Prestwich et al. (2013 7"’ Hopkins et al. (2002), Sage et al. 

(1992], 

Table 1. Distances, oxygen abundances, SFRs, X ray luminosities, and observed number of HMXBs for the galaxies in our sample. 
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pend strongly on L±. Hence, we took a large value so that for 
most galaxies L t h < L,h, 2 . and the number of sources IV 38 above 
Lth ,2 is determined by simple counting. Only for 26 of our 49 
galaxies is L t h > L t h ,2 and then a correction based on the XLF is 
needed, 


N 38 =N + q SFR L“ dL. (2) 

'-'Ah,2 

For 22 of them, the integral in the right-hand side of Eqn. [2] is 
<s 1, and N 38 = N. In the few remaining cases, the correc¬ 
tions amount to 33% (UGC 5720), 50% (SBS 0335-052), 86% 
(NGC 4194), and 100% (IIZw 40) of the final value. The lumi- 

G 38 G 36 

nosity L x was determined in the same way as L x , but using 
Zth ,2 instead of L lhii in Eqn. IT] The corrections in this case are 
lower than 7%, except for NGC 4194 (15%). 

The corrections made by Eqns. [1] and [2] allow us to obtain a 
uniform sample of the numbers of sources and luminosities of 
the galaxies. It is important to note that most of these corrections 
are small, and that they do not affect our final results, as N 38 , 

G 38 G 36 

L x , and L x ' are used only to visualize and derive preliminary 
trends. The full Bayesian analysis of the data, from which our 
final results are obtained and our conclusions are drawn, relies 
on the original data (N, L x , and L th ). 
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Fig. 1. Oxygen abundances calculated from raw spectral data by using 
both fits to the N 2 ratio given by |Pettini & Pagel] ( |2004[ >, as a function of 
the abundances taken directly from the literature. Cyan dots (red stars) 
correspond to the linear (cubic) calibration of the N 2 index; the dashed 
line represents the identity. 


2.2. Star formation rates 

We extracted the star formation rates (SFRs) of 19 of the galaxies 
of our sample from Mineo et al.| ( |2012| . These authors use a 
combined estimator for the SFR based on the near-UV (NUV) 
radiation from young stars, corrected by the IR emission due to 
heating of dust, 


SFR = SFR 


-NUV 


■(1-i/)SFRir, 


(3) 


where q is a correction factor related to the IR emission com¬ 
ing from the old stellar population, whose value is null for 
starburst galaxies and 0.4 for normal disk galaxies. The NUV 
and IR contributions to SFR (in M 0 yr 1 ) can be calculated 
from the 2312A and 8 - 1000 pm luminosities (in ergs -1 ) as 
SFR ir = 4.6 x 10 44 T, r and SFR NUV = 1.2 x 10 -43 L NUV , re¬ 
spectively. In addition to the galaxies taken from the work of 
|Mineo et al.| ( |2012[ ), we adopted this estimator for those galaxies 
for which NUV and IR observations are available. For the rest 
of the galaxies we used far-UV (FUV) SFRpuv estimations cor¬ 
rected by a factor of 1.23. According to Prestwich et al. (|2013|> 


Brorby et al. (2014), this factor gives the ratio of FUV esti 


and 

mations to those based on NUV and IR observations. For a single 
galaxy (II Zw 40) the only available SFR estimation is based on 
H ff measurements. In this case, we transformed this value into a 
FUV SFR using the conversion of Hunter et al. j( 2010| >, and then 
once again using the factor 1.23 to transform it to the NUV/IR 
SFR scale. In this way, we obtained a sample as homogeneous 
as possible to compare the effects of SFR in the production of 
HMXBs. 


2.3. Metallicities 

The metallicity estimations for our sample were taken directly 
from the literature. As we are dealing with star-forming galax¬ 
ies, in which emission lines are used to estimate abundances, we 
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use the oxygen abundances 12 + log(0/H) as a metallicity indi¬ 
cator. The abundance values compiled are a priori far from ho¬ 
mogeneous, as they were obtained with different methods whose 


systematic differences are still unclear (Bresolin et al. 2004 

Gar- 

nett et al. 2004; Liang et al. 2006 

Yin et al. 2007am Kew 

ey & 

Ellison 2008; Lopez-Sanchez et a 

. 2012). Although it is still not 


certain which method is better to estimate the absolute oxygen 
abundance of a galaxy, it is believed that the relative metallicities 
between galaxies estimated with the same calibration are reliable 
( jKewley & Ellison|2008[ >. However, due to the large metallicity 
range of the galaxies in our sample, it was not possible to calcu¬ 
late the oxygen abundances with a single method because of the 
lack of homogeneity of the spectral lines observed. 

It is important to note that in this work we are not interested 
in the computation of absolute oxygen abundances, but in the 
correlation of metallicity with the X ray emission of galaxies. 
With this aim and in order to evaluate the impact of comparing 
abundances estimated by different calibrations, we searched the 
literature for the raw spectral information of each galaxy, and 
built a homogeneous subsample of abundances computed by a 
single method. To this aim we used the intensities of the spectral 
lines [Nil] A 6583 and H„ , and computed the abundances of 21 
galaxies using the N 2 index ( |Pettini & Pagel|2004| ). This method 
is single-valued, and is based on a ratio between lines that are 
very near in the spectrum, hence it does not depend on the red¬ 
dening correction or flux calibration. 


For those 21 galaxies of our homogeneous subsample, we 
calculated the abundances both using the linear and cubic fits 
for the N 2 ratio given by Pettini & Pagel| ( |2004| >. The correlation 
between both N 2 estimations and the abundances taken directly 
from the literature are shown in Fig. [T] The variance between 
them is 0.17 dex in the case of the linear method and 0.22 dex 
for the cubic fit. Hence, despite the possible systematic devi¬ 
ations between the different calibrations, the general trend ob¬ 
served supports the use of the complete sample to estimate rela¬ 
tive metallicities. 
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Fig. 2. X-ray luminosity as a function of SFR for the sample of galaxies 
compiled in this work. Black circles correspond to those galaxies with 
measured X-ray luminosities, while blue triangles represent the detec¬ 
tion limits for single sources. These can be regarded as upper limits for 
the galaxy luminosity in the case it had only one source. The red dashed 
line (log L x = log SFR + 39.4) is the relation found by |Mineo et al.| 
( [20121 ). It describes well the data at high SFRs, whereas it fails to do so 
in the opposite SFR range. For log SFR < —1.1, all the data fall above 
the fit, which suggests the existence of metallicity effects in the sample. 


3. Analysis of the sample 

3. 1. Observed trends 


In Fig. [2] we show the Lx-SFR relation for the galaxies in 
our sample, together with the best fit of Mineo et al.| (2012 [ 
log Lx = log SFR+39.4). In order to make a meaningful compar¬ 
ison, in this plot we use L^’ 36 so that all galaxies have the same 
luminosity threshold, which is also the one used by these au¬ 
thors. It is evident that their best fit describes properly the mean 
behavior of high-SFR galaxies, mostly taken from their work, 
but fails to describe the correlation at low SFRs. In this limit, 
the X-ray luminosities are higher than those predicted by the fit. 
Indeed, all galaxies with log SFR < —1.1 fall above the fit of 
|Mineo et al.| ( |2012| ). The fit also fails to describe the large dis¬ 
persion of the data in the whole metallicity range. As the SFR 
of galaxies is correlated to their mean metallicity, the increase in 
the X-ray luminosity with respect to the fit of Mineo et al. ( 2012[ > 
observed at low SFRs could be due to metallicity effects. 

However, it could also be due to statistical fluctuations, as 
low-SFR galaxies have small HMXB populations. The increase 
in the dispersion at low SFRs suggests that this effect is present 
in our sample. Disentangling both effects is crucial to reach 
a meaningful conclusion about the metallicity dependence of 
HMXB populations. 

To further investigate the influence of metallicity on the Lx- 
SFR relation, in Fig. 3 we show the dependence of the L^' 36 /SFR 
ratio on the oxygen abundance of the galaxies. Despite the high 
dispersion, the observed behavior is consistent with an anticor¬ 
relation between the X-ray luminosity per unit SFR and the oxy¬ 
gen abundance. It is also apparent that a linear relation is a good 
description of the correlation. The best linear fit (not including 
upper limits) is logL^ 36 = 39.26 + log SFR - 1.01 [log(0/H) - 
log(0/H)o l, taking 12 + log(O/H) 0 = 8.69 from As£lund_et 


al. 


2009). The L^' 36 /SFR ratio is similar to that of 


Mineo et 


Fig. 3. X-ray luminosity per unit SFR as a function of the oxygen abun¬ 
dance of the galaxies in our sample. Black filled squares correspond to 
those galaxies with measured X-ray luminosities, while blue triangles 
represent the detection limits for single sources. These can be regarded 
as upper limits for the galaxy luminosity in the case it had only one 
source. The red short-dashed line represents the value found by |Mineo| 


et al. 


(2012) for solar-metallicity galaxies (L°’ 36 /SFR = 39.4), while 


r G,36 


the gray dotted line is the best linear fit to our sample, logL^ - 
39.26 + log SFR - 1.01[log(O/H) - log(0/H) 0 ]. A clear trend of de¬ 
creasing L^' 36 /SFR with oxygen abundance is observable, which is a 
direct evidence of the existence of metallicity effects in HMXB popula¬ 
tions. Stochasticity effects are also evident, and require a more refined 
treatment to make a robust prediction of the variation of the L^' 36 /SFR 
ratio with metallicity. The cyan and green long-dashed lines above and 
below 12 + log(0/H) = 8 respectively, represent our best models for the 
high- and low-metallicity subsamples, with the 68% confidence bands 
for the parameter Q (small-dotted gray lines). 


al. (2012 [ ) at solar metallicity, while it increases about 1 dex per 
dex of decrease in abundance. However, the dispersion around 
the regression line is still important (~ 0.5 dex), suggesting that 
stochasticity effects should be investigated to make a more robust 
prediction of the variation with metallicity of the X-ray luminos¬ 
ity per unit SFR of galaxies. 

The increase in the L X /SFR value toward low oxygen abun¬ 
dances might be due to an increase either in the HMXB pop¬ 
ulation size, or in the mean HMXB luminosity. To disentangle 
these two physical effects, in Fig. [4] we show the observed num¬ 
ber of sources A 38 with luminosities above the fixed threshold 
Lth ,2 = 10 38 erg s 1 , per unit SFR, as a function of the metallic¬ 
ity of the galaxies. A clear trend is seen, in the sense that this 
number is lower for high-metallicity galaxies (log(0/H) > 8) 
than for low-metallicity ones (log(0/H) < 8). The difference is 
about 1 dex. The presence of many data points with error bars 
much smaller than this difference suggests that it cannot be ex¬ 
plained by statistical effects. As the observed number of sources 
depends mainly on the population size, and only weakly in the 
HMXB luminosities, this result suggests that the main metallic¬ 
ity effect is a change in the total number of HMXBs per unit 
SFR. Indeed, the observed difference in A 38 /SFR from low- 
to high-metallicity galaxies is roughly the same amount as the 
change in L^ ,36 /SFR derived from Fig. 3 This implies that the 
increase in the integrated luminosity of galaxies at low metallic- 
ities can be explained by the increase in the size of their HMXB 
populations. An increase in the mean luminosity of HMXBs, if 
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Fig. 4. Observed number of HMXBs with luminosities higher than 
10 38 erg s _1 , per unit SFR, as a function of the oxygen abundance. Black 
circles are the data points, with error bars representing Poissonian un¬ 
certainties. Large error bars reaching the lower end of the plot represent 
galaxies with no detected sources; their upper ends correspond to A 38 = 
1. Our best model for low-nretallicity galaxies (12 + log(0/H) < 8, 
green solid line, Q = 3.13, cr extIa = 0) predicts ten times more HMXBs 
than the corresponding best model for high-metallicity galaxies (cyan 
dot-dashed line, Q = 2.14, er cxtra = 22). For both models, the 68% con¬ 
fidence bands for the parameter Q and the dispersion predicted for the 
data are shown (gray dashed lines and shaded regions, respectively). 
The mean Q = 2.28 for high-metallicity galaxies with cr extra = 0 is 
also plotted (blue solid line). The difference between low- and high- 
metallicity galaxies cannot be explained by statistical effects. 


any, is certainly much lower. Once again, the dispersion of the 
/V 3X /SLR data is high, indicating the presence of stochasticity 
effects which must be taken into account to obtain a more robust 
prediction of the amount of the observed variations with metal- 
licity. 


3.2. Stochasticity effects 

As pointed out before, all the data shown in this section present 
large dispersions. Part of them may be due to the stochasticity 
of the observed HMXB populations, which affects both their 
sizes and observed luminosities. Indeed, it has been shown by 
Gilfanov et al. ( |2004j > that the collective luminosity of a small 
population of randomly chosen sources can grow non-linearly 
with the number of sources, contrary to the prediction of models 
based on infinite populations. To obtain a robust description of 
the trends present in the data, we developed a simple stochastic 
model for the number of sources and total luminosity of the pop¬ 
ulation of HMXBs of a star-forming galaxy with a given SFR 
and oxygen abundance, which attempts to describe the observed 
behavior (see Appendix |A| for its details). Following Brorby et 
|al. | P014| l, we assume that the mean number of sources in a 
galaxy follows a Poisson distribution whose mean scales with its 
SFR, 


(Nhmxb) = 6 SFR, (4) 

where g is a free parameter that may depend on the oxygen 
abundance. We also assume for the HMXBs a standard power- 
law XLF with index a = -1.6 in the range [L mm , L max ] = 
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[10 35 ,10 40 ] erg s" 1 ( |Fabbiano|2006l|Mineo et al.|2012| . The ob¬ 
served number of sources and integrated luminosity of a par¬ 
ticular galaxy above any specified threshold luminosity can be 
computed using the XLF. We note that our model includes only 
the dependence of the population size on metallicity, and not 
that of HMXB luminosities, as the XLF is fixed. In addition, 
our model includes two fluctuation sources, an intrinsic one due 
to the stochasticity of the star formation process, and another 
one given by the detection process, which selects only those bi¬ 
naries with luminosities above L t h- Model fitting was done by 
a Markov-chain Monte Carlo method, which gives the best-fit 
value of g, together with its variance. 

( 2014[ ) proposed that g is different for 


As Brorby et al. 


low- and high-metallicity galaxies, we divide our sample into 
two subsamples comprising galaxies below and above 12 + 
log(0/H) = 8, respectively. We fit our model to each sample, 
obtaining log glow = 3.13 ± 0.13 and log ghigh = 2.28 + 0.02 
(which correspond to % w = 12.8 and g^gh = 1-81 using the 
notation of |Brorby et al. |20 1 4|p' This results agree reasonably 
well with those of |Mineo et~ak ( 201 2[ > (q = 1.49) and [Brorby] 
et al. ( |2014 1 (q = 14.5), and imply that low-metallicity galax¬ 
ies produce ~ 10 times more HMXBs per unit SFR that their 
high-metallicity counterparts. 


The dispersion given by our model describes the one of the 
low-metallicity subsample very well, as shown in Fig. [4] How¬ 
ever, this is not true for the high-metallicity data. To investi¬ 
gate this issue, we added an extra dispersion cr extra to the num¬ 
ber of sources predicted by our model, and fit it to the high- 
metallicity subsample for different (fixed) values of cr ex t ra . A 
Bayesian model comparison shows that the data are best de¬ 
scribed by cr extra = 22, typically 3 times that given by Poissonian 
effects, and a slightly lower log ghigh = 2.14. The decrease in the 
latter value compared to the <x extra = 0 model is due to the fact 
that the extra dispersion (constant and larger than the Poissonian 
uncertainties of the data) increases the weight of lower-A' values 
relative to higher-A' ones, hence lowering the best-fit value. The 
fact that an extra dispersion is not observed in the low-metallicity 
sample could be explained by the different sample sizes (556 
high- vs. 12 low-metallicity sources). Indeed, the test of picking 
random samples of 12 high-metallicity sources and performing 
the same analysis shows that the extra dispersion remains unde¬ 
tected if the size of the high-metallicity sample is small. More¬ 
over, the mean best-fit value of ghigh obtained for the random 
samples of 12 high-metallicity sources agrees to within uncer- 
taintes with the best-fit value of the whole high-metallicity sam¬ 
ple. This fact implies that the increase in the number of HMXBs 
per unit SFR at low metallicities is not an effect of the difference 
in the sample sizes. 


The aforementioned results confirm that the metallicity de¬ 
pendence observed in our sample is not an effect of statistical 
fluctuations, but a real physical trend, as the value of g is dif¬ 
ferent for the low- and high-metallicity subsamples even when 
the dispersion of the data is well described. Indeed, the probabil¬ 
ity that the low-metallicity data comes from random fluctuations 
around the best high-metallicity model (including the extra dis¬ 
persion) is ~ 10 47 times that given by the best low-metallicity 
model. Hence, in agreement with Brorby et al. ( |2014 i, our anal¬ 
ysis suggests that the number of HMXBs per unit SFR increases 
by a factor of ~ 10 at low metallicities, even if we take into ac¬ 
count the statistical fluctuations. As an interesting by-product. 


1 The factor Q in our work differs from that used by Brorby et al. 
1 20141 in the normalization of the luminosity function. 


















































V. M. Douna et al.: Metallicity dependence of HMXB populations 


these fluctuations are well described by Poissonian effects at low 
metallicities but are larger at high metallicities. 

The models that best fit the number of sources of the galaxies 
also describe nicely their integrated luminosities (Fig. [3]». In the 
high metallicity limit, the predictions are also in well agreement 
with the relation of Mineo et al. ( 2012j >. This result suggests that 
the most numerous HMXB populations produced at low metal¬ 
licities can explain not only the different observed number of 
sources, but also the enhanced integrated X-ray luminosity of 
galaxies. Hence, the dependence of the instrinsic luminosity of 
the sources with metallicity would be small. However, given that 
the observed number of HMXBs depends on the XLF through 
the detectability of the sources, and that the XLF is fixed in our 
model, a dependence of the luminosity on metallicity cannot be 
completely discarded. 

To test the consistency of the observed luminosities with 
those predicted by our model in a statistical way, we extended 
our model to predict the luminosities of the observed HMXBs 
by sampling the XLF with the restriction L > L± (see Ap¬ 
pendix [A] for details). We fit the model to the high- and low- 
metallicity subsamples, now including the observed luminosi¬ 
ties of the galaxies together with their numbers of HMXBs in 
the comparison. We included the previously obtained extra dis¬ 
persion in the high-metallicity model. The resulting values of the 
free parameter Q are in full agreement with the ones obtained be¬ 
fore, confirming that the luminosity dependence on metallicity is 
small, if any. 

To further explore this possibility, we analyzed the correla¬ 
tion between the observed mean luminosity of HMXBs in each 
galaxy (L^' 38 /A 38 ) and its metallicity (Fig. 5b. This figure shows 
a very weak trend (red, dashed line) in the sense that HMXBs 
are brighter in low-metallicity populations, masked by large- 
amplitude fluctuations. If this trend could be confirmed by more 
numerous and more precise data, then it might indicate that the 
observed variation of the total X-ray luminosity of the HMXB 
populations results from a large contibution due to the change in 
the population size, plus a minor one originated in the variation 
of the mean HMXB luminosity. 

Finally, it is important to point out that although the simplest 
model that describes the dependence of the number of HMXBs 
on the oxygen abundance is a split function, there is no physi¬ 
cal evidence to justify the break at 12 + log(0/H) - 8. Hence, 
following the apparent linear trend shown by the number of ob¬ 
served HMXBs in Fig. [4] as a second approach we also param¬ 
eterized the factor Q by a linear function of the oxygen abun¬ 
dance. A fit to our whole sample confirms the existence of a 
dependence of Q on metallicity with a confidence greater than 
1 - 5 x 10 \ However, the best-fit parameters suggest that the 
anticorrelation between the number of binaries and the metallic¬ 
ity of the galaxies is milder than that predicted by a sharp break. 
In this case the change in the population size per unit SFR is 
only -0.14dex per dex. The discrepancy is originated in the bet¬ 
ter statistics of high-metallicity data, which have more weight in 
the fit than low-metallicity ones, hence biasing the slope to shal¬ 
lower values. A deeper exploration of this issue would require 
larger samples, and a more detailed treatment of population size 
fluctuations. 


4. Conclusions 

We compiled from the literature a large sample of 49 galax¬ 
ies with measured star-formation rates, oxygen abundances and 
HMXB population properties (observed number of sources and 
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8.5 
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Fig. 5. Mean luminosity of HMXBs above a threshold of 10 38 ergs -1 , 
as a function of the oxygen abundance. Black, filled squares are the 
data points, while the red dashed line is the best linear fit to them, 
log/IV 38 = (38.7 ± 1.3) - (0.18 ± 0.16)[log(O/H) - log(0/H) 0 ], 
A low-significance anticorrelation between these variables is observed, 
which might suggest a brightening of HMXBs at low metallicities. As 
in previous figures, the dispersion around the fit is large (0.4 dex). 


luminosities above some threshold luminosity). This sample en¬ 
hances those previously compiled by Min eo et al.| ( |2012| ) and 
|Brorby et al. | ( |2014| ) by giving the oxygen abundances of the 
galaxies, and by including both the number of sources in each 
galaxy together with its integrated luminosity. As the low- and 
high-metallicity galaxies were taken from different works, with 
different uncertainties and selection effects, the data were ho¬ 
mogenized to make statistical comparisons and draw meaning¬ 
ful conclusions. The SFRs in our sample were corrected to a 
common estimator. For the oxygen abundances, we used a sin¬ 
gle estimator in the largest possible subsample and the rest of 
the sample was shown to be unbiased with respect to this homo¬ 
geneous set. Any residual non-uniformity that may have been 
left over by this process is small, and does not affect our conclu¬ 
sions. Hence, our complete sample is then suitable to investigate 
the effects of metallicity in the properties of HMXB populations. 

Our main results are the following: 


The size of the HMXB populations of low-metallicity (12 + 
log(0/H) < 8) galaxies, per unit SFR, is ~ 10 times larger 
than that of solar metallicity galaxies. This result is robust; 
it describes well the observed number of HMXBs and in¬ 
tegrated X-ray luminosities of nearby star-forming galaxies, 
taking into account the observed fluctuations. It does not de¬ 
pend on the different sizes of the low-and high-metallicity 
subsamples. It is also in excellent agreement with those of 
Kaaret et al. ( |2011| and|Brorby et al. ( 2014[ ), who found 


similar results using slightly different methods and data. 

The dispersion in the observed number of HMXBs of low- 
metallicity, low-SFR galaxies can be explained by Poisso¬ 
nian fluctuations associated with the initial mass function 
and XLF sampling. This is not the case for high-metallicity 
galaxies, which show fluctuations typically half an order of 
magnitude higher than those expected from these effects.[Mi]] 
neo et al. (2012 ) > have reported an analog dispersion in the 


Lx-SFR relation of the high-metallicity subsample, which 
cannot be attributed to selection effects. Our results indicate 
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that the large dispersion is not an effect of the IMF or XLF 
sampling. Therefore it could be due to a different physical 
origin, most probably the internal metallicity dispersion of 
the galaxies and the age distribution of HMXBs. Indeed, nu¬ 
merical simulations of stellar populations which compute the 
effect of these fluctuations give dispersion values that agree 
with our findings (Artale et al., in prep.). This effect is hidden 
in the low-metallicity subsample behind the large Poissonian 
fluctuations that arise from the small number of sources. 

- Present evidence suggests that the aforementioned metal¬ 
licity effects occur below 12 + logfO/H) « 8. Although a 
smooth linear trend in the correlations is apparent, it is not 
supported by statistical tests on the present sample. A sharp 
break at this metallicity provides the best description of our 
data. However, the smoothness of the metallicity dependence 
could be affected by differences in the sample size or residual 
inhomogeneities. A larger sample including galaxies near the 
break and more low-metallicity galaxies is crucial to solve 
this issue. A more thorough description of the data fluctua¬ 
tions, especially those from non-Poissonian sources, would 
also be desirable. 

- The dependence of the luminosity of HMXBs with metallic¬ 
ity, if any, is much smaller than that of the population size. 
However, refined models are needed to determine its signifi¬ 
cance. 


Our results support the general picture describing the en¬ 
hancement of the HMXB population size at low metallicities 


(Dray 2005 Linden et al. 2010 

Kaaret et al.|2011| 

Basu-Zych| 

et al. 2013a|b| Fragos et al.|[2013 

Brorby et al.~p014 

). This im- 


plies that these energetic sources were more numerous in the 
Early Universe, hence a promising source of heating and partial 


reionization of the IGM, as proposed by several authors (Power 

et al. 

2009; Mirabel et al.||2011| Jeon et al][2013j (Power et al. 

2013 

Knevitt et ai.|2014). It also implies that the energy feed- 


back from HMXBs to the ISM could play an important role in 
the formation and evolution of low-mass galaxies, supporting re¬ 


cent results ( [Artale et al.|2015j l. To assess the exact magnitude of 
these effects, a more detailed description of the variation of the 
number of HMXBs and their luminosities must be given. This is 
only possible with further deep observations that unveil the low- 
luminosity end of the XLF in nearby galaxies, together with ho¬ 
mogeneous metallicity and SFR estimations. Age measurements 
of the star-forming regions would also be valuable to investigate 
the effects of this important parameter. 
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Appendix A: Monte Carlo model for HMXB 
populations 

Appendix A. 1: The model 


Our model aims at describing the observed number of sources 
and total luminosity of the population of HMXBs of a star¬ 
forming galaxy with a given SFR and oxygen abundance, to 
perform a robust comparison with observations. We start from 
the fact that the progenitors of HMXBs are massive stars, and 
that these sources have short evolutionary times (less than ~ 
lOOMyr; ^Belczynski et aLl|2004] [20081 |20T0l [Shtykovskiy &| 
Gilfanov 2005 200 7|l. Hen ce, their number should correlate with 
the SFR, as Brorby et al. ( |2014j ) propose in their Eqn. 4. How¬ 
ever, star formation can be regarded as a stochastic process in 
which the final result (the formation of a single or binary star, 
and in the last case the initial masses of the components, the 
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semimajor axes and eccentricities, etc.) follows a certain proba¬ 
bility distribution. The sampling of this distribution in a partic¬ 
ular stellar population is a source of fluctuations in the number 
of HMXBs (Ahmxb) of a galaxy. If we assume that the proba¬ 
bility Phmxb of producing a HMXB is fixed, then Ahmxb has 
a Bernoulli distribution. As the total number of stars A* in any 
stellar population is large, and Phmxb ^ 1 because the initial 
mass function is a strongly decreasing function, we can approx¬ 
imate Ahmxb as a Poisson variable. The mean value of Ahmxb 
is then (Ahmxb) = A*phmxb- However, as it is not our aim to 
compute the distribution of Ahmxb from first principles, follow¬ 
ing Brorby et al. (20141 we parameterize it as 


(Ahmxb) = G0og(O/H)) SFR, 


(A.l) 


as A* scales with SFR. The factor Q is the same as q in the Eqn. 4 
of|Brorby et al. | (2014 1, apart from the normalization factor of 


the XLF. We note that the difference between our model and that 
of these authors is that the above equation decribes the mean 
number of sources in a Poisson process, while the corresponding 
equation of Brorby et al. ( |2014| describes the exact number. 
When N ~ 1 (as in low-SFR galaxies) differences arise, and 
the stochasticity must be taken into account. Our factor Q may 
depend continuously on the oxygen abundance, or just assumed 
to be different for low and high metallicities. 

A second source of stochasticity is given by the detection 
process, which selects only those binaries with luminosities 
above some threshold L ltl . The XFF sampling in the HMXB 
population introduces statistical fluctuations in the number of 
sources detected. Given that detection is a Bernoulli process, the 
composition with the Poisson process described above produces 
a new Poisson process with mean 


<A de t> = Pdet0(log(O/A)) SFR, (A.2) 

where 


Pdet ““ 


t'l i 


i//(L)dL. 


(A.3) 


Here the XFF is — k norm L for L e [T lnm , T max ], and k norm — 
[‘ . L a dL is its normalization factor. 

-'f^min 

Our model allows, for a fixed set of free parameters, the com¬ 
putation of the number of sources detected in a galaxy Adet, by 
sampling a Poissonian with mean given by Eqn. A.2 The lu¬ 
minosities of these sources are computed numerically by sam¬ 
pling the XFF Adet times, with the restriction that L > L±. The 
sum of the sampled luminosities is Ldet , the total luminosity of 
the sources detected in the galaxy. In this way, we can com¬ 
pute a realization of the observables (Adet, Aiet) for any galaxy 
with known SFR and 12 + log(0/H). We stress that our model is 
stochastic, i.e., each realization of our model gives a different set 
(Adet, Tdet)- Therefore, the comparison with observations must be 
made with care. We also stress that it is possible for our model 
to predict that no sources will be detectable in a given galaxy 
(Adet = Tdet = 0), hence upper limits in the observed luminosity 
of galaxies can be included in the comparison. 


to observations and fix the values of the free parameters; instead, 
we use a Bayesian approach. We note that the free parameter of 
the model may be Q itself, or any set introduced to parameter¬ 
ize Q as a function of metallicity (e.g., the slope and intercept 
of a linear function). We will refer to them collectivelly as 6. 
From Bayes’ theorem, the posterior probability density / post of 
the model parameters 6, given the observed data ({A, L^}, with 
i e [1,49]) is 


/ post (e|{A,Z°}) 


£({N,L°}\6)f pn (e) 
f £({N,L°mf pti (0)d0’ 


(A.4) 


where f pr i is the prior probability density of the parameters, and 
£ the likelihood of the data given the values of the parameters. 
The best-fit parameters are those that maximize the posterior 
probability density / post . 

Assuming that we have no previous information on the pa¬ 
rameters, we adopt a uniform prior. On the other hand, the like¬ 
lihood is the product of the individual joint probability distribu¬ 
tions of the number and luminosities of HMXBs of each galaxy. 


49 

-£({A, z°})|0) = Y] 6)P N (Ni\6), (A.5) 

i=i 

where fi is the conditional probability density for the observed 
luminosity given the observed number of sources, P N the proba¬ 
bility of the observed number of sources, and the product arises 
because the data of different galaxies are independent. As de¬ 
scribed in the previous section, Te best fit gives in this case the 
values 


P N (Ni\6) = 


e A^><A det )M 

A,-! 


(A.6) 


Assuming that the measurement uncertainty for the luminosity 
L ( x '. of galaxy i is Gaussian with standard deviation cr L i , 


fdPi.m = 


vw L , 


r 

Jo 


exp 


^dct) 

2cr h 


, 2 \ 


,gL(Tdet|A/)r/Tdet- 

(A.7) 


Here gzXTdetlA) is the probability density of the luminosity of 
the detected HMXBs. This factor has no analytical expression. 


but the integral in Eqn. A.7 can be computed numerically by a 


Monte Carlo method as the mean of the Gaussian in the inte¬ 
grand over a large number M » 1 of realizations of our model, 
i.e.. 


M 


- 


y/2ncr L jM 


X exp 


(^X,i ^obs ,j) 

2cr 2 r ■ 

L,i 




(A.8) 


In this way, the likelihood £ and hence the posterior probability 
density / post can be computed within our model. The maximiza¬ 
tion of the latter was performed by a Markov-chain Monte Carlo 
scheme, giving the best-fit parameters and their uncertainties. 


Appendix A.2: Model fitting 

The stochasticity of our model prevents us to use classical statis¬ 
tical methods such as least-squares fitting to compare the model 


Article number, page 9 of [ 9 ] 




























